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We discuss a novel damping mechanism of r-mode oscillations in neutron stars due to processes 
that change the number of protons, neutrons and electrons. Deviations from equilibrium between the 
number densities of the various species lead to the appearance in the Euler equations of the system 
of a dissipative force, the so-called rocket term. Such a force affects the evolution of the r-mode 
oscillations of a rotating neutron star and we estimate the corresponding damping timescale. In the 
description of the system we employ a two-fluid model, with one fluid consisting of all the charged 
components locked together by the electromagnetic interaction, while the second fluid consists of 
superfluid neutrons. Both components can oscillate and we consider two different kind of r-mode 
oscillations, one predominantly associated with comoving displacements of the two fluids and a 
second one associated with countermoving, out of phase, displacements. In the former case we find 
that the dissipation mechanism associated with the rocket term is not effective in preventing the 
growth of r-mode instabilities. In the latter case the rocket term prevents the growth of the r-mode 
instability for temperatures larger than about 10 9 K. In our analysis we include the mutual friction 
dissipative process between the neutron superfluid and the charged component. In order to simplify 
the computation we neglect the interaction between the two r-mode oscillations as well as effects 

Oh' related with the crust of the star. Moreover, we use a simplified model of neutron star assuming a 

y> , uniform mass distribution. 

PACS numbers: 04.40.Dg, 04.30.Db, 26.60.-c, 97.10.Sj, 97.60.Jd 
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I. INTRODUCTION 



> 

o 

Rapidly rotating neutron stars have been the subject of intensive investigations in the last years. Of particular 

interest are neutron star oscillations, which might be useful to shed light on the internal structure of these stars [l|, 0). 

t ' There are various modes of oscillations, among them r-mode oscillations are probably the most interesting ones, 

because they provide a severe limitation on the star's rotation frequency through coupling to gravitational radiation 

^ ■ emission 0, [j] . These oscillations can be damped by various dissipative process that can take place in the star |5|, |(| . 

However, if dissipative phenomena are not strong enough, these oscillations will grow exponentially fast in time until 

the star slows down to a rotation frequency where some dissipative mechanism sets in. Since neutron stars are observed 

to rotate at very high frequencies, any model of neutron star must provide an efficient mechanism of dissipation of 

r-mode oscillations. In this way, the study of r-mode damping is useful in constraining the stellar structure and can 

d ' be used to rule out some exotic phases of matter [7| . 

Standard neutron stars are stellar objects with a mass of about 1.4M and a radius of about 10 Km. They are 
believed to have a crust of about 1 Km, with an outer part made of a lattice of ions embedded in a liquid of electrons 
and an inner part made of nuclei embedded in a liquid of 1 So superfluid neutrons. In the interior of the star nuclei 
are melted and both neutrons and protons are expected to condense into BCS-likc supcrfluids. However, neutron 
interaction in the 1 Sg state at supernuclear matter density is repulsive, but it is still possible to form Cooper pairs 
in the 3 P2 channel [8(. The proton density is much smaller than the neutron density and they can form pp Cooper 
pairs in the isotropic Sq channel. Pairing between protons and neutrons does not take place because of the large 
mismatch between their Fermi energies. In the core of neutron stars muons might be present (when fi e > m^), or 
deconfined quarks in a color superconducting phase, moreover pion or kaon condensates might be realized. In the 
present paper we shall not consider any of these possibilities and assume the core of the neutron star comprises only 
neutrons, protons and electrons. 

R-modc oscillations have been studied extensively in the literature and various damping mechanisms have been 
proposed [l| . Shear and bulk viscosity are able to suppress the instability in two different ranges of temperatures 0, @ . 
For temperatures smaller than about 10 5 K shear viscosity will always suppress the instability, but with increasing 
temperature the effect of shear viscosity is gradually suppressed. For very high temperatures, bulk viscosity becomes 
an efficient mechanism for damping r-mode oscillations. In the superfluid phase bulk viscosity is suppressed by Pauli 



blocking, but at temperatures above than 10 10 K -10 11 K standard nuclear matter is believed to be in the normal 
phase and for temperatures larger than 10 10 K bulk viscosity is able to damp the instability. Therefore, one has 
an "instability window" for standard neutron stars corresponding to a range of temperature approximately given by 
10 6 K - 10 10 K. The exact values depend on the details of the model considered. The instability window is in part 
reduced by the "surface rubbing" between the core and the crust of the star jjj [U| ■ This mechanism results in a 
viscous boundary layer between the core and the crust of the star which damps r-mode oscillations for temperatures 
less than about 10 10 K and for sufficiently small frequencies. 

In Rcfs. [Il| - [l4| it is studied the effect of mutual friction in reducing the instability window. It is shown that the 
typical timescale of mutual friction is of the order of 10 4 s and is therefore too long for damping the r-mode instability. 
Indeed, the timescale associated with gravitational wave emission is of the order of few seconds (for a millisecond 
pulsar). However, mutual friction can reduce the instability for certain values of the entrainment parameter |12| or 
for large values of the drag parameter [14| . 

In the present paper we study a novel dissipative processes associated with the change in the number of protons, 
neutrons and electrons. In real neutron stars these processes can take place in the outer core and in the inner crust of 
the star and are related to beta decays and interactions between the neutron fluid and the crust. Both these processes 
lead to the appearance of a dissipative force (the so-called rocket term) in the Euler equations of the system. This 
force is due to the fact that when two or more fluids move with different velocities a change of one component into 
the other results in a momentum transfer between the fluids. This change in momentum is not reversible, because 
it is always the faster moving fluid that will lose momentum. The resulting dissipative force is proportional to the 
mass rate change, and to the relative velocity between the fluids. The name "rocket term" reminds that the same 
phenomenon takes place in the dynamical evolution of a rocket whose mass is changing in time as it consumes its 
fuel. As far as we know, the dissipative force due to the rocket term has not been considered in the context of r-mode 
oscillations. In the mass conservation laws, it is in general assumed that the neutron, proton and electron numbers 
are separately conserved quantities. 

In order to simplify the analysis we consider a simplified model of neutron star consisting of a fluid of neutrons, 
protons and electrons and no crust. Protons and electrons are locked together by the electromagnetic interaction and 
therefore we consider that the system consists of two fluids. As a further simplification we assume that the star has a 
uniform mass distribution with density p = 2.5po and a radius of 10 km. Since this simplified model of star does not 
comprise a crust we consider only number changing processes associated with weak interactions. 

In our analysis we consider two different r-mode oscillations. One is associated predominantly with toroidal comoving 
displacements and the second one is dominated by toroidal countermoving displacements. We shall refer to these 
oscillations as "standard" r-modes and as "superfluid" r-modes, respectively. These two modes decouple for a star 
made by uniform and incompressible matter, and we shall restrict to treat such a case. We find that the dissipative 
mechanism associated with the rocket term is not efficient in damping standard r-mode oscillations. On the other 
hand, this mechanism becomes efficient in damping the superfluid r-mode instability at temperatures of the order of 
I0 9 K. Therefore this mechanism becomes efficient before bulk viscosity starts to be able to damp r-mode oscillations. 
In agreement with [1 31 ] we find that the superfluid r-mode oscillations are damped as well by the mutual friction force 
for sufficiently large values of the entrainment parameter. 

This paper is organized as follows. In Scction[H]we review the hydrodynamic equations in the two fluid approxima- 
tion. In the Euler equations we neglect shear and bulk viscosities but consider the effects due to the mutual friction 
force and to the rocket term force. In Section IIIII we present the differential equations governing the deviations from 
equilibrium of a rotating two fluid system. We study the linearized problem, neglecting deformations of the star due 
to rotation and assuming that the star has a uniform density. Although not very realistic, these assumptions allow 
us to simplify the study of the r-modes, while concentrating on the strength of the effect of the rocket term. The 
contribution of Urea processes to particle mass creation rates are evaluated in Section IIVI In Section [V] we evaluate 
the timescale corresponding to the rocket term force and compare the result with the timescale associated with other 
dissipative processes. We draw our conclusions in Section IVII In Appendix [A] we report some details about the 
evolution equations for the comoving and countermoving displacements. The expressions of the decay widths of Urea 
processes in presence of various nn condensates are reported in Appendix [Bj 

II. HYDRODYNAMICAL EQUATIONS FOR THE MULTIFLUID NEUTRON STAR MODEL 

The equations describing the dissipative processes of neutrons, protons and electrons in the outer core of a stan- 
dard neutron star have been studied in detail in Rcf. |ll| . In general, the entropy production rate depends on 19 
independent coefficients which are related to the various dissipative processes. However, for vanishing temperature 
and neglecting viscosities one obtains the expression given in Ref. [l5[ which depends only on two different coeffi- 
cients. One is related to mutual friction and the second one with the so-called rocket term. Here we review the basic 



hydrodynamical equations in presence of these two different dissipative mechanisms, neglecting the presence of shear 
and bulk viscosities. 

For a system consisting of neutrons, protons and electrons, the mass conservation law is given by, see e.g. [111. Il5j. 

dtp x + Vi{p=ovi) = r x , (i) 

where T x is the particle mass creation rate per unit volume and the index x = n,p,e refer to the particle species, that 
is, neutrons, protons and electrons. In these equations we have considered that some process can convert neutrons 
in protons and electrons and vice versa. Therefore, we are assuming that the various components are not separately 
conserved. One possible mechanism leading to a change in the particle number densities is given by the weak processes 

n — >• p + e~ + v e , p + e~ — >• n + i/ e . (2) 

These reactions lead to a change in the chemical potentials of the various species and therefore are associated with 
number density changes. 

A different process is the so-called "transfusion", see Ref. [16[. In this process when a compression takes place, the 
ionic constituent of the crust are squeezed and part of their neutron content is released and augments the supcrfluid 
neutron component. The opposite mechanism, related to a reduction of the pressure leads to the neutron capture by 
the ions of the crust. In the present section we consider that a generic mechanism is at work to produce a change 
in the number densities. In Section IIVI we shall evaluate the particle mass creation rate corresponding to the beta 
decay processes. Regarding the transfusion processes the corresponding creation rates are difficult to evaluate and we 
postpone their calculations to future work. 

In any case, the three particle creation rates are not independent quantities, because charge conservation implies 
that 



whereas baryon number conservation leads to 



= r P , (3) 



r P = -r„ , (4) 



meaning that only one creation rate is independent. 

It is possible to simplify the treatment of the system considering that our analysis regards processes that happen at 
a time scale much larger than the electromagnetic time scale. Therefore, we can consider that electrons and protons 
are locked together to move with the same velocity [13], see however [18 |. Moreover, charge neutrality implies that 
the number densities of electrons and protons are equal, i.e. n e = n p . Therefore, electrons and protons can be treated 
as a single charge- neutral fluid and henceforth we shall refer to this fluid as the "charged" component, employing the 
subscript c to label it. As a matter of fact, the system can be viewed as consisting of two fluids, with mass densities 

pn = m n n„ and p c = m„ n c , (5) 

where m n = m p + m e and n e = n p = n c . The Eulcr equations obeyed by the two fluids are given by 

fMF 

(dt+vfrjXv? + e„w i ) +V i (/x„ + $) + eW i ^ = ii— , (6) 

Pn 
fMF p 

(d t + viV 3 ){vi - e cWl ) + ViO«c + $) - e<VV 4 ^ = - J -±— + (1 - e„ - e c ) — w % , (7) 

J Pc Pc 

where i,j label the space components, we have defined a chemical potential by mass p x = p x /m n , and w = v c — v„ 
represents the relative velocity between the two fluids. The quantities e n and e c are the entrainment parameters, that 
are related to the fact that momenta and velocities of quasiparticles may not be aligned [19| . 
The gravitational potential, <&, obeys the Poisson equation 

V 2 $ = 47rG(p„ + Pc ) . (8) 

The last term on the right hand side of Eq. ([7]) is the rocket term. This force is due to the fact that when two 
fluids move with different velocities a change of one component into the other results in a variation of the momentum 
of each fluid component. This change in momentum can be viewed as a force proportional to the mass rate change, 
r„, and to the relative velocity between the two fluids, w. Actually, in Eqs. (|6l7p . one can see that the rocket term 
acts only on the charged component. The reason is that part of the rocket term contribution can be included in the 



definition of the mutual friction force; here we have employed the same notation used in Rcf. [T^]. In the analysis of 
the possible dissipative mechanisms of star oscillations this term is usually neglected. Indeed, it is in general assumed 
that the neutron, proton and electron numbers are separately conserved quantities, that is T p = T e = T n = 0. 

The force term /j MF entering into both Euler equations corresponds to the mutual friction force between the two 
fluids. At the microscopic level it is due to the scattering of electrons off the magnetic field trapped in the superfluid 
vortices [lj}. Indeed the core of the superfluid vortices is magnetized, because of the interaction between the neutron 
and the proton condensates. The mutual friction force can be expressed as 

jf F = 2p n B'e ijk il j w k + 2p n Be ijk n j e klm fl l w m , (9) 

where the coefficients B, B' can be written as 

72 1Z 2 

B = -, and B' = = , (10) 

1 + 1Z 2 ' 1 + ll 2 ' v ' 

where 1Z is the dimensionless "drag" parameter |ll| . The value of 72 depends on the mechanism of trapping of the 
magnetic field inside superfluid vortices [171 120( . The details of this mechanism are unknown and one can consider 
three different regimes: the weak drag regime, 72 -C 1, the strong drag regime 72 3> 1 and the intermediate drag 
regime, 72 ~ 1. For small values of the drag parameter one can express the coefficients B, B' as a function of the 
cntrainment parameter. According with Ref. [ll]], one has that 

B = 4 x 1(T 4 ^L= ( ^-)" ( £ A and B' ~ B 2 . (11) 

vT^lo.05/ \10 14 g/cm 3 y 

III. PERTURBED HYDRODYNAMICAL EQUATIONS 

A non-vanishing mass creation rate influences the evolution of the various hydrodynamical. Indeed, the continuity 
equations ((T|) as well as the Euler equations (|7|) depend on r„. Therefore, in the analysis of the various modes of 
oscillations of a neutron star one has to take into account effects related to this term. In the present paper we only 
discuss its effect on the evolution of the r-modes of a superfluid neutron star, although it would be equally interesting 
to study its effect on other pulsation modes. In the following analysis of the hydrodynamical equations we also include 
the mutual friction force and we follow very closely the recent a naly sis of the r-mode oscillations developed in Ref. [21j 
for normal fluid stars, and extended to superfluid stars in Ref. [14|. 

As in Rcf. [14J . we study the linearized hydrodynamical equations for the perturbations around an equilibrium 
configuration of a neutron star rotating with constant angular velocity fl, the r-modc is purely axial to leading order 
in n and we work in the equilibrium hydrodynamical quantities, and we assume that the background configuration 
is such that the two fluids move with the same velocity, thus at equilibrium w = 0. 

It is useful to write the Euler equations for the perturbed quantities using as degrees of freedom (dof) the center 
of mass displacement and the relative displacements between the neutron fluid and the charged fluid. We define the 
comoving velocity as 

Sv = ^5v n + ^<5v c , (12) 

P P 

and the countermoving, or relative, velocity as 

(5w = (5v c — <5v n . (13) 

The continuity equation for the comoving degree of freedom is not affected by the rocket term and is given by 

a t 8p + V j (p5v j ) = 0, (14) 

on the other hand the continuity equation for the countermoving dof depends on it. It is useful to employ as a second 
continuity equation the one for the charged fraction, x c = Pc/ p, which is given by 

i r„ 

dtSx c = — V ■ [x c {\ — x c )p(5w] — Sv ■ Vx c . (15) 



The rocket term affects the Euler equations for both the comoving and countermoving velocities and one has that 

d t Svi + 2e ijk n j 6v k + -V Z 6 P - %V lP + V;<5$ = (1 - e) — S Wi , (16) 

P P 2 P 

d t (l - e)6wi + \7 l {5(3) + 2B'e ljk n j Sw k - 2Be ljk ^e klm riiSw m = (1 - e)—Sw l , (17) 

where here we have defined I = e c + e„ = e„(l + p n / Pc) = ^n/x c , and where 

5(3 = 5fi c - 5fl n (18) 

and 

B = B/x c , B' = 1 - B'/xc . (19) 

The hydrodynamical equations can be studied employing a perturbative expansion of the various hydrodynamical 
variables in H, the star rotation frequency. Actually, the expansion is in the parameter Q/Qk, where fix is the Kepler 
frequancy of the star. For superfluid systems, this expansion is particularly convenient, as one can show that the 
complicated system of equations for the comoving and countermoving degrees of freedom decouple as these variables 
depend on different powers of Q. 

For our study we consider some simplifying, admittedly unrealistic, assumptions. We neglect the deformation 
of the star due to rotation, which affect the different hydrodynamical variables at order O 2 . We use the Cowling 
approximation, that is, we neglect perturbations of the gravitational potential associated with the oscillations of 
the star. As a further simplification, we also consider a model where the mass density of the star is uniform. As 
emphasized in the Introduction, our goal is to study the effect and impact of the rocket term in the evolution of the 
r-modes, and we leave for future studies a more realistic model of the star. 

Oscillations of a fluid clement of a stars can be described by the Lagrangian displacement vector £, which can be 
decomposed into a sum of toroidal and spheroidal components. Since neutron stars can be described employing the 
two fluid model, one defines comoving and countermoving displacements, respectively £+ and £_, by means of the 
equations 

5v = <9 t £+cxO£ + , 8w = d t $- cxft£-. (20) 

These two displacements describe the center of mass oscillation and the out of phase oscillation of the two fluids, 
respectively We then expand these quantities in terms of toroidal and spheroidal components 

£+ = r S~] I 0, -r^dtf,, ~I<i m d 9 ) Y im +rV Si m ,Zi m d 8 , -r^A ) Y lm , (21) 

l,m l.m 

i- = r yV 0, -r^Tjd'P, -hmde j Yi m + r y^ [si m , zi m dp, -t^A j Yim , (22) 

l,m l,m 

where Yi m are the spherical harmonics. The fluctuations of the pressure and of the chemical potential difference can 
be written respectively as 

$P = Pgr^2(i m Y lm , (23) 

l.m 

6/3 = grY,r lm Y lm , (24) 

I .11! 

where g = Vt^r (with fig = GM/R 3 ) fixes the scale of pressure and chemical potential fluctuations. Notice that with 
these definitions, i\ m and Cim are dimensionless. 

Since in a superfluid star one has two different kind of displacements, in principle one can have two different kind of 
r-mode oscillations, one associated with the comoving dof and one associated with the countermoving dof. Actually, 
the hydrodynamical variables defined above obey a complicated set of coupled differential equations, see |14| . with 
couplings between comoving and countermoving displacements. However, as shown in [14|, at the leading order in ft, 
one finds that the equation for the comoving displacement decouples and one can determine an analytic expression 
for the standard r-mode oscillation. Regarding the mode associated with the countermoving dof, it turns out to be a 
general inertial mode. That is, it is not a mode dominated by the toroidal components. However, for incompressible 
stars with uniform density one has that this inertial mode turns into an r-mode. We shall restrict to this case and 
analyze this r-mode oscillation in Sec. IIIIB1 We explicitly consider the effect of the mutual friction in the equations 
of motion, the reason is that in this way we can analyze the regime where the mutual friction coefficients, B and B' , 
are large. Therefore our results will explicitly depend on the values of these parameters. 



A. Standard r-mode oscillations 

For the standard r-mode oscillations one assumes that the comoving toroidal displacement , Ki m , is of order unity, 
while the spheroidal comoving displacements arc of order £7 2 . All the countcrmoving displacements turn out to be 
of order f2 2 as well. Since the standard r-mode oscillation is dominated by Ki m , it is very similar to the r-mode 
oscillation in normal fluids [l| , and can be easily determined after imposing proper boundary conditions [141 ] . To first 
order in the rotation frequency of the star, one has that the typical frequency of the oscillations (measured in the 
corotating frame) is 



2mfi 

1(1 + 1) 



(25) 



In our analysis we restrict to analyze the case I = m = 2, which corresponds to the most unstable r-modc. 

Regarding the pressure perturbations, they are of order f2 2 , whereas 8(3 oc H 4 [12, LLJ]. The order in Q of the toroidal 
oscillations and of the pressure and chemical potential fluctuations are reported in the first line of Table I. 
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0(Q 2 ) 


o(n 2 ) 

O(Q ) 
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0(Q 2 ) 


o(n 2 ) 

o(n 2 ) 



TABLE I: Order in SI of the comoving and countermoving displacements, of the pressure fluctuation and of the chemical 
potential fluctuation for the standard r-mode oscillation and for the superfluid r-mode oscillation. 



For purposes of estimating the damping time scales associated to both mutual friction and the rocket term that we 
will carry out in Sec. [Vj wc have to determine the solutions for the countermoving dof. The equations governing the 
evolution of the various dynamical variables are reported in the Appendix [X] Assuming constant mass density and 
hydrostatic equilibrium we find that t; +1 obeys the following differential equation 



2 // 

r T i+i 



(A x +B t - l)rr/ +1 + (A 2 B 2 - A 1 B 1 ) n+1 - A 2 B 4 



J+i 



R 2 



(26) 



where the prime indicates differentiation with respect to r and the coefficients Ai , Bi are reported in Appendix [X] As 
shown there, the last term on the right hand side of this equation arises because we have assumed that matter is in 
hydrostatic equilibrium. The differential equation has solution given by 



T l+1 (r) = f(r) + C x r n ^+W 



(27) 



where /(r) is the particular solution of the differential equation and where C\ and C 2 are the coefficients of the 
homogeneous solution, to be fixed by the boundary conditions. The exponents of the homogeneous solution are given 
by 



"■1,2 



Ai+Bi± yffAi + B x y + 4(A 2 B 2 - A X B X ) 



(28) 



and it turns out that n 2 is negative, meaning that in order to avoid divergences at r = 0, it must be C 2 = 0. It is 
interesting to note that for vanishing mutual friction one has that n\ =1 — 1 and n<i = — (l + A). As a second boundary 
condition we assume that the chemical potential difference vanishes at the surface of the star, that is S/3(R) = 0. 

For completenes we report the equation obeyed by the radial component of the countermoving spheroidal displace- 
ment, which is given by 



c = 



r/ +1 AxrTi +1 



A, 



Ao 



(29) 



B. Superfluid r-mode oscillations 



Assuming that k\ m is of order unity one finds that the spheroidal countermoving displacements are of order f2 2 . The 
driving force on the countcrmoving displacement is the chemical potential difference which turns out to be of order 
fl 2 . The order of the toroidal comoving displacement depends on the compressibility of the fluid. For a compressible 



fluid it is of order il°, while for an incompressible fluid it is of order f2 2 . The reason can be traced back to the fact 
that comoving oscillations are driven by pressure oscillations and it turns out that Ki m ~ f2 _2 For compressible 
fluids the pressure oscillations are proportional to chemical potential oscillations and therefore £ ~ Vt 2 and thus Ki m 
must be of order unity. Moreover, for this kind of mode, comoving spheroidal displacements turn out to be of the 
same order in tt of comoving toroidal displacements, meaning that for a compressible fluid this oscillation is a generic 
inertial mode and not an r-mode. Since for a compressible fluid various components of the displacements are of the 
same order in f2, one has to solve a system of coupled differential equations. 

The situation is much easily tractable for incompressible fluids. In this case one can assume that spheroidal 
oscillations arc of order 0(il 2 ) and then toroidal oscillations turn out to be of the same order. The order in $7 of the 
various displacements and of the pressure and chemical potential fluctuations for incompressible matter are reported 
in Table I. We shall restrict the analysis to the case of incompressible fluids, where the comoving and countermoving 
dof decouple, with the superfluid r-mode oscillation dominated by the toroidal displacement k\ m . To first order in the 
rotation frequency of the star and to first order in the entrainment parameter, the typical frequency of the superfluid 
r-mode oscillation (measured in the corotating frame) is 

a :d + e). (30) 



' 1(1 + 1)' 

As for the standard r-mode, we restrict to analyze the case I = m = 2, which corresponds to the most unstable 
r-mode. Moreover we consider only small values of the entrainment. 

The analysis of the the superfluid r-mode oscillation is very similar to the one we have performed for the standard 
r-mode oscillations, with the roles of Ki m and of the pressure oscillations interchanged with ki m and the chemical 
potential oscillations. We find that for superfluid r-modes, ki m obeys the same equation that Ki m obeys for standard 
r-modes, and the chemical potential fluctuation obeys the same equation that pressure fluctuation obey for standard 
r-modes. Regarding the pressure oscillation Q m one has to solve an equation analogous to Eq. (|26p . but without the 
last term on the right hand side, because we are now considering an incompressible fluid. We find that 

0+i = C ir Sl + C 2 r s * , (31) 

where s\^ depend on the parameters of the model. One of the two coefficients is always negative, and therefore in 
order to avoid the divergence at the origin, we have that 

O+i - Cr s . (32) 

We fix C by demanding that the comoving toroidal displacement, Ki m , is properly normalized, as in Ref. [13| . 

IV. URCA PROCESSES AND THE ROCKET TERM 

In order to evaluate the damping timescale associated with the rocket term we need to evaluate the mass creation 
rate T n . There are two main processes that are responsible of the change of the densities of neutrons, protons and 
electrons in neutron stars. The capture/release of neutrons from the ions in the crust and the Urea processes. The 
first process regards the properties of ions at extremely high densities and it is difficult to take into account properly. 
One should consider all the microscopic processes as well as macroscopic processes (like crust-quakes) that can lead to 
a transfusion of material between the crust and the underlying fluid. We are not aware of any calculation of the mass 
creation rate associated with the transfusion mechanisms. However, it is worth mentioning that such processes might 
become important for a self-consistent treatment of r-modes. R-mode oscillations have a radial component that can 
be of the order of 100 m. Such a radial displacement comes with a pressure oscillations that can lead to a change in 
composition of the crust and may lead to transfer of material between the crust and the underlying superfluid. 

We postpone the treatment of the transfusion processes to future work. For the time being we restrict to Urea 
processes that take place among protons, neutrons and electrons. It was found in Ref.[22[, that for certain realistic 
equations of state the direct Urea processes are allowed when the star density exceeds the nuclear saturation density 
po = 2.8 x 10 14 g cm -3 , and the proton fraction exceeds the threshold value x p = \. We shall restrict to consider 
such processes assuming that these conditions are realized in the interior of neutron stars. Henceforth we shall take 
x p = Xp , but our results do not strongly depend on the precise value of x p as far as x p > x c . The direct Urea reactions 
are given by 

n — » p + e~ + v e , p + e~ — >• n + v e , (33) 



and the corresponding reaction rates are given by 

d 3 p 



r£rca=/ It ^/-(l^plfl-JelEMW^-i 5 /)- (34) 

i—n,p,e,iy L ^ ' - 1 spins 



i— n,p,e,u 

for the neutron decay process, and by 

r& rca = / II \ (2k)22e] fpfe (1 ~ fn) E \M\ 2 (27r) i 6^(P i - P t ) , (35) 

i— n,p,e,f spins 

for the electron capture process. In both reaction rates, \M\ 2 is the squared of the scattering cross section for the 
weak interaction process (see Refs. |22| and [23| for its explicit value), and 

f* = e (Bi-MO/r - 1 ' (36) 

is the Fermi-Dirac distribution function of the particle specie i = n,p, e. 

In equilibrium, or quasi-equilibrium, the neutron decay and the electron capture processes balance, meaning that 
furca = -^urca- I n sucn a case j one nas r n = T p = 0, and the rocket terms in the hydrodynamical equations (|T6|) 
and (]17p vanish. At equilibrium the velocities of the proton and neutron fluids should be the same, and further, the 
chemical potentials should fulfill the relation 

fin = Me = Mp + Me ■ (37) 

Therefore, a perturbation around the equilibrium distribution is related to a difference in the velocities of the two 
fluids and/or a deviation from chemical equilibrium. Out of equilibrium the neutron decay process and the electron 
capture process do not compensate, resulting in a net source (or sink) of neutrons, and thus, in a non-vanishing value 
of T„, given by 

r„ = r Urca — r Urca . (38) 

In order to evaluate this quantity we consider small perturbations of the equilibrium distribution functions, fi — fi+Sf; 
upon substituting the perturbed distribution function in Eqs. (|34l35j) and considering terms linear in the perturbation, 
we obtain the mass creation rate 

r„ = -— (S[i n - 5fi c + -^-(1 - e« - e P )(<5w) 2 J f Urca , (39) 

in agreement with the results of Ref. [15| . Note that the last term in Eq. (|39|) gives the correction to the chemical 
equilibrium due to the relative motion between the two fluids. 

The value of Turca depends on the particular neutron supcrfluid phase that is realized in the star. In Ref. [23j three 
phases have been considered, one corresponding to a Sq condensate and two corresponding to the ^Pi condensate with 
uij = and rrij = 2. In any superfluid phase the Urca reaction rate is suppressed with respect to the corresponding 
value in unpaired matter. The reason is that the phase space available for scattering is reduced in the superfluid 
phase. This leads to a suppression of the neutrino emission rate and of the bulk viscosity coefficient in the cores of 
superfluid neutron stars [23|. In our case this leads to a suppression of the rocket term force. However, it is worth 
remarking that while the suppression is exponentially large for the So and for the 3 P2 condensate with rrij = (of 
order exp(— A/T), with A the corresponding energy gap), it is much smaller for the 3 P2 condensate with rrij = 2. 
This is due to the fact that in the latter case the quasiparticle fermionic dispersion law is gaplcss. We report the 
results of Ref. [23[ regarding these three superfluid phases in Appendix [Bj 

In the situation where the direct Urca processes are not allowed, one should consider modified Urca processes, 

n + N ->• p + N + e~ + v e , p + N + e~ ->• n + N + v e , (40) 

where N is an additional nucleon. However, the reaction rates of the modified Urca processes are 3-5 orders of 
magnitude lower than those of direct Urca processes [24| , and the effect of the rocket term in the damping of r- modes 
will be weaker. For the time being, we will assume that the conditions for the direct Urca processes are allowed in 
the star. 

In general, one can express the mass creation rate as 

58 w 
r„ = f^— S, (41) 

k B T ' y ' 



where ks is the Boltzmann constant and we have defined 

S/3™ =m n 6j3-^{l- e„ - e c )(5w) 2 , (42) 

where in the case when the Urea processes are allowed S = Furca and its value depends on the superfluid phase 
considered. 

V. STABILITY ANALYSIS 

Given the mass creation rate, we can now evaluate the damping timescale associated with the rocket term and 
compare the result with the timescale of different processes. In principle one can determine the damping times 
associated with any mechanism from the imaginary part of the mode frequency. However, in the present analysis we 
shall limit ourselves to the energy integral estimates. 

First we evaluate the energy associated with any oscillation of the star. In the frame corotating with the star, see 
[13| | for more details, it can be written as 

|H 2 + (1 - e)x p (l - x p )\Sw\ 2 ) + ]T VijtfHdpl - ^^ ] dV , (43) 

i,j=n,c J 

where the integral is extended to the volume of the superfluid core of the star and we have defined V%j = 4^ . The 
energy stored in the oscillations contains a kinetic energy contribution, Ek , and a potential energy contribution, E pot . 
The variation of the total energy due to mutual friction and rocket term dissipations can be obtained using the 
continuity and Euler equations for the oscillations, integrating by parts, and discarding surface terms. In this way we 
have that the total energy variation associated with the oscillations is given by 

BE _d(E k +E pot ) _ ,_.,., , , , ™ 2 - 




j dV\\Swf M F\ + \Sn 2 ^) ■ (44) 

These results are in agreement with the general considerations for the entropy generation due to both mutual friction 
and the rocket term of Ref. (15j . Further, the requirements of positive entropy production allow us to guarantee that 
these terms always represent energy losses, and not gains. 

The damping timescale associated to the mutual friction is then given by 

i=-2i(fL=^/* |iw - f "'' (45) 

while for the rocket term one has that 

1 1 fdE\ 1 



^ 2E\6tU 2Ej dV ^^- ™ 

The timescale associated with gravitational wave emission is instead given by 

7- = -^(lP) =^£^f +1 (|AJ 2 + |JU 2 ), (47) 



2E\dtJ gw 2E 
where uj r is the r-mode frequency and w, = lo t — m Q is the frequency measured by an inertial observer. Moreover 



4ttG (Z + !)(/ + 2) 
' 1(1 -l)[(2l + l)\l] 



n = ^ ,,, K ,^;::L* , m 



and Di m and Ji m are the mass and current multipoles respectively. 

The three timescales tmf, trt and r gw depend on the particular r-mode oscillation considered. In the following 
sections we shall separately analyze standard r-modes and superfluid r- modes. In order to make our results comparable 
with the analysis of Ref. [12J we define the entrainment parameter 

(49) 



l—x c — e 

and according with [12| we consider the range of values e < 0.06 
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FIG. 1: Damping time tr T associated with the rocket term as a function of the entrainment parameter e (defined in Eg, (1491) ) 
for the weak drag regime. We have taken T = 10 9 K and a critical superfluid temperature T c — 10 10 K. The star has uniform 
density p — 2.5po, with radius R — 10 km and mass M ~ 1.47M©. 



A. Timescales for standard r-modes 

The standard r-mode oscillation is predominantly associated with toroidal comoving displacements. According with 
Table I and equations (|20p we have that Sv oc ft, while Sw oc f2 3 to leading order in the Q expansion. Thus the leading 
term in the total energy, Eq. (|43p . is associated with the comoving oscillations and therefore E oc f2 2 . 

We can now evaluate the power counting in f2 of the various timescales. For the mutual friction timescale, we have 
from Eq.© that /mf ~ £l 4 and it follows from Eq. (|4"5l) that tmf oc £l~ 5 . As for the damping time of the rocket term 
in Eq. (|46]) . we have that trt oc $7 -6 . This comes from the fact that 5(5 oc J7 4 at the leading order in 51. This also 



means that in Eq. (|4"2"| one can neglect the second term on the right hand side, because (<5w) 2 oc f2 6 . The timescale 
associated with gravitational- wave emission for standard r-modes has been computed in Refs. [J, y, 0, M, IIH, l26| ■ For 
a star of constant density, the growth timescale of the I = 2 current multipolc is given by 



'gw 



22 



M 



lAMp 







R 



lOkr 



P 

lms 



and therefore r gw oc fl~ 6 . 

In r-mode analysis it is useful to introduce the frequency independent quantities t mf , 
from the corresponding timescales extracting the Q dependence. For standard r-modes we have that 



.0 



RT' 'gw 



(50) 



r" , which are obtained 



1 

tmf 



1 



'MF 



n 



\J~xGp 



1 
trt 



1 

'rt 



O 



y/nGp 



1 



1 



fi 



\/irGp 



(51) 



Since we have explicitly considered the mutual friction force in the Euler equations, the mutual friction timescale, 
t mf , and the rocket term timescale, t^ t , depend on the drag parameter. In agreement with Ref. [14( we find that for 
R ~ 1 the mutual friction is able to damp the r-mode oscillations, while for large or small values of R the damping 
time is too long. 

Regarding the rocket term damping timescale, we compute the numerical value of the coefficient 3 in Eq. (|4ip . In 
Fig. [T] we report the plot of the rocket term damping timescale versus entrainment. We consider the weak drag regime, 
TZ <C 1, and that nn pairing takes place in the 3 P2{mj = 2) channel, which has the smallest reduction factor (reported 
in Appendix [B] and corresponding to case C). However, results with different values of 1Z and with pairing in the 
1 5'o channel and in the 3 P2(m,j = 0) lead to very similar results. We assume throughout that the critical temperature 
of the superfluid phase is T c = 10 10 K. We find that r^ T does weakly depend on the entrainment parameter. The 
damping is stronger for small values of the entrainment, but in any case the damping time is extremely large. 

More in detail, the rocket term damps the r-mode instability for 



trt < T g 



(52) 



and considering that the trt has the same frequency dependence of T gw , it follows that the condition for damping is 
independent of the rotation frequency of the star and can be written as 



'RT 



<r„ u 



(53) 
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Thus, when this condition is realized, the star is rotationally stable for any value of f2, as far as the Q/Qk expansion 
makes sense. 

However, r° is of order of tens of second, and as is clear from Fig. Q] the timescale of the rocket term is way 
too big for damping the r-mode instability. More precisely, considering the temperature dependence of the rocket 
term, we find that for temperatures below 10 15 K the damping timescale of the rocket term is larger than the growth 
timescale of gravitational waves. Therefore the effect of the rocket term is completely negligible in the analysis of 
standard r-mode instabilities. The qualitative reason why the rocket term effect is negligible is related to the fact 
that it is associated with countermoving displacements, which are subleading in standard r-mode oscillations. The 
suppression in f2 can only be compensated by a very efficient particle conversion mechanism. In the present analysis 
we have assumed that particle conversion is due to direct Urea process, which is the most efficient microscopic decay 
one can consider. Therefore, the rocket term can have a significant effect on standard r-mode oscillations only if a 
very efficient macroscopic (transfusion) mechanism of particle conversion is at work in the neutron star, weak decays. 

B. Timescales for superfluid r- modes 

The superfluid r-mode oscillation for incompressible supcrfluids is dominated by toroidal countermoving displace- 
ments. According with Table I, we have that Sw ex f2 and Sv (x ft 3 , to leading order in the Q expansion. The total 
energy in Eq. f|43[) is now dominated by the countermoving displacement and the power counting in Q is the same 
one has for standard r-modes, i.e. E ex f2 2 . 

Regarding the mutual friction timescale, we have from Eq. (|9]) that /mf ~ Q 2 ■ Therefore in the case of superfluid 
r-modes the mutual friction force is much larger than in standard r-modes. This is clearly due to the fact that mutual 
friction force is proportional to the countermoving displacement, which are dominant for superfluid r-modes. From 
Eq. (|45[) we have that tmf oc f2 _1 , which is four orders in Q smaller than the corresponding expression for standard 
r-modes. 

As for the damping time of the rocket term, we now have that 5/3 ex fi 2 and therefore from Eq. (j4l))) we have that 
trt oc fl~ 2 . Also in this case the timescale is much shorter than for standard r-modes, and for the same reason. 

The timescale associated with gravitational-wave emission for superfluid r-modes is much larger than the corre- 
sponding value for standard r-modes. The reason is that gravitational wave emission is only associated with comoving 
displacement which are suppressed with respect to the countermoving modes. As a matter of fact the mass and 
current multipolcs turn out to be both of order Q 3 and we find that r gw ex Q~ w . Thus, for superfluid r-modes we 
have that 



i i / n \ i i / n \ 2 i i / n 



tmf t mf \^ttGpJ ' t rt r^ T \y/TiGpJ r gw r° \y/nGp 



(54) 



We shall now analyze in detail the stability of the superfluid r-modes. In addition to mutual friction and rocket term 
one has to consider shear and bulk viscosity as well. We have evaluated both of these quantities for superfluid r-modes 
and found that 



1 1 / fi \ 1 1 / n 

t sv t s ° v \^ttGpJ ' T bv r° v WnGp 



(55) 
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T gw t sv Tbv tmf trt 





In the numerical evaluation of the bulk viscosity timescale we have employed the results of |23[ , where the bulk 
viscosity is evaluated for direct Urea processes, but taking into account the reduction factor associated with the 
presence of the superfluid phase. 

The study of the stability of r-mode oscillations now requires to take into account all the different timescales 
simultaneously. The relation for the stability is given by 

i i 

(56) 

meaning that when this relation is satisfied one has a critical condition for stability. When the quantity on the left 
hand side is negative, then the mode is unstable. 

In Fig. [2] we report the result for the superfluid r-mode "instability window" for a star with uniform density, 
p = 2.5/9o, and R = 10 km. The dashed red line represents the instability window in the absence of the rocket term 
and mutual friction. The region above the dashed red line is unstable with respect to gravitational wave emission 
when only shear and bulk viscosity damping mechanisms are considered. At low temperature, shear viscosity is the 
dominant dissipative mechanism. With increasing temperature shear viscosity is less efficient and it can damp r-mode 
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FIG. 2: Instability window of the superfiuid r- modes of a star with uniform density p — 2.5po, with radius R = 10 km and mass 
M ~ 1A7Mq. In both panels, the dashed red line represents the instability window in the absence of the rocket term. The 
full blue line represents the instability window with the inclusion of the rocket term. The left panel corresponds to the case of 
vanishing entrainment and therefore vanishing mutual friction. In the right panel we consider two values of entrainment. The 
horizontal full lines correspond to the effect of the mutual friction for e = 0.0002, lower line, and e = 0.002 upper line. In our 
simplified model of star, the mutual friction is independent of the temperature. 



oscillations for smaller and smaller values of the frequency. The behavior of bulk viscosity is the opposite and starts to 
damp r-mode oscillations for temperatures of the order of 10 10 K. The full blue line represents the effect of the rocket 
term. It has a behavior qualitatively similar to the bulk viscosity, but it becomes effective at smaller temperatures. 
Therefore, the instability window with the inclusion of the rocket term is much reduced. In the left panel of Fig. 
[2] we have assumed zero entrainment and hence, vanishing mutual friction. In the right hand panel of Fig. [2] wc 
include the effect of mutual friction. In the simplified model of neutron star that we are considering, mutual friction 
is independent of the temperature. We consider only the case of small drag parameter, where B and B' are given in 
Eq. tfTTj) and are functions of the entrainment. The lower horizontal line corresponds to e = 0.0002, while the upper 
line corresponds to e = 0.002. In agreement with [l3|, we find that for values of e ~ 0.02 the effect of the mutual 
friction is to damp superfiuid r-modes for any value of the frequency. 



VI. CONCLUSION 



Superfiuid neutron stars are characterized by various oscillation modes. Of particular interest are r-mode oscilla- 
tions, because they are unstable via gravitational wave emission, meaning that they would lead to the spin-down of 
the star, unless some dissipative mechanism sets in and damps these oscillations. 

We have considered a very simple model of star comprising a fluid of neutrons, protons and electron, with uniform 
density and no crust. We have assumed that protons and electrons arc locked together by the electromagnetic 
interaction and therefore the system can be described as made by two fluids: the neutron superfiuid and the charged 
component. 

In this model one has to consider two different r-mode oscillations. The standard r-mode oscillation, which are 
predominantly toroidal comoving displacements of the two fluids, and the superfiuid r-mode oscillations which are 
associated with toroidal countermoving displacements. In realistic neutron stars these two modes are coupled, however, 
in the limit of small rotation frequency and assuming that the star has a uniform density and is incompressible, they 
decouple. 

We have evaluated the effect of the rocket term on the stability of these two r-mode oscillations. In our analysis 
we have included the effect of mutual friction and of viscosities and compared the corresponding timcscales with the 
timescale of gravitational wave emission. We find that the dissipative mechanism associated with the rocket term 
is not efficient in damping standard r-mode oscillations. The reason is that standard r-modes are predominantly 
comoving modes, while the rocket term is associated with the countermoving flow, which is subleading. On the 
other hand, this mechanism becomes efficient in damping the superfiuid r-mode instability. The reason is that in 
this case the countermoving mode is the leading oscillation of the system and we find that for temperatures of the 
order of 10 9 K the rocket term has a typical timescale comparable with the timescale of gravitational wave emission. 
Quite interestingly, for superfiuid r-modes the rocket term becomes an efficient damping mechanism before bulk 
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viscosity sets in. In agreement with [13[ we find that the superfiuid r-modc oscillations can be damped as well by the 
mutual friction force for sufficiently large values of the entrainment parameter. Unlike the mutual friction force, the 
rocket term is not strongly dependent on the entrainment parameter and it reduces the instability window even for 
vanishing entrainment. Since the entrainment parameter is a poorly known coefficient, the rocket term represents an 
interestingly mechanism for superfiuid r-mode damping. 

It would be interesting to extend the present analysis to more realistic equation of states, in order to see whether 
the effect persists. It would be also quite interesting to study the effect of "transfusion" of neutrons and protons from 
the crust to the outer core. For standard r-modes we find that the damping timescale of the rocket term is quite large, 
however this timescale would be reduced if the transfusion of nucleons from the crust to the core results in a mass 
rate larger than the one due to weak processes. In this case the rocket term might be able to reduce the instability 
window for standard r-modes as well. 

We have limited our analysis to the region of small chemical potential fluctuations, 5(3 <T. As in Ref. [27j , it would 
be interesting to extend the present analysis to the suprathermal region, 8/3 > T, in order to study the spind-down 
time of unstable stars. The rocket term might also be relevant for hybrid stars with a quark matter core. In that 
case one should consider weak processes between quarks of different flavors and transfusion processes between nuclear 
matter and quark matter. 

Appendix A: Evolution equations 

Assuming uniform density the comoving and countermoving modes decouple |14j . For a star with uniform density 
one can impose hydrostatic equilibrium obtaining 

P(r) = G^(R 2 -r 2 )p\ (Al) 

where we have assumed that the pressure vanishes at the surface of the star. Here R = 10 km is the radius of the 
star and we shall consider a mass density p = 2.5po, where po is the saturation density of nuclear matter. With these 
values we obtain that the mass of the star is M ~ 1.47M Q , where Mq is the mass of the sun. In this case we have 
that the pressure and the various components of the countermoving mode obey the following set of equations |14j 

km = asi + i+bz l+1 (A2) 

t 1+1 = ch m + dzi + i + esi + i (A3) 



Or, 



i+i 



Or 



= -2t ;+ i - / ki m - g zi+i + h sj+i (A4) 



-3s l+1 - — Q +1 + p zi+i , (A5) 



dr ^ T 1 

where the various coefficients have been derived in Ref. |14| and for I = m are given by 



B-iB' 1 

a = - : = ^ ; (A6) 

(1 - e) - B> - iB V2T+3 

b = (l + 2)a (A7) 

(, +2 ) (1 - e -)-m'- iS l±±M±f±^ ( A9) 



; + 2 V 5 + 2/ 



^ 



. IB' - iB- — l -] (A10) 

1 + 2 V 5 + 21 J y J 

f ■■ -(l+l)c (All) 

g = -pe (A12) 



h = uj z 1(1 - e) - ^ B -^jfj (A13) 

p = (l + l)(l + 2), (A14) 



where Co = Si/Op and w = a/flo, with f2o = \/GM/R 3 . In Eq. (|A5[) we have that 



V= 9 -^ and r=^, (A15) 

P dlogp ' v ; 
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which depend on the equation of state. Since we assume hydrostatic equilibrium we have from Eq. (|A1I) that 

2r 2 

V = -^ tt and T = 2 . (A16) 

R z — r z 

Finally, according with Ref. jl4| , we have that the pressure fluctuations are given by 

^ = 2 vm—iR^- (A17) 



Upon substituting the expressions above in the equations (|A2|) . (|A3|) . ()A4[) and (|A5j) and expressing ki m and zi+\ in 
terms of s/ + i and t; + i we have two coupled differential equations for si+\ and 77+1. These equations can be written 
as a second order differential equation 



„2_" 



{A 1 + B x - l)rr/ +1 + (A 2 B 2 - A 1 B 1 )t 1+1 - A 2 B 4 -, j , (A18) 



where 



ad ~ be ac + e , . » 

51 = - 3 -^' (A21) 

52 = 5Td' ( A22 ) 

i? 4 = r -^— . (A23) 

The analysis of the superfluid r-modes is analogous to the one we have done for the standard r-modes. However, in 
order to have an r-mode oscillation and not a generic inertial mode, one has to assume that the fluid is incompressible 
PJJ . In this case r — > 00 and the differential equations one has to solve are simpler. 

Appendix B: Reaction rates for the Urea processes in different superfluid phases 

The reaction rates for the direct Urea processes that may occur in superfluid neutron stars depend on the particular 
pp and nn superfluid condensates that can be realized. Various superfluid phases have been considered in Ref. |23| . 
and the corresponding reaction rates have been evaluated; here we review some of their results. 

If there is no superfluidity, the reaction rate for the direct Urea process can be expressed as 

f Urca - (A/)1.667 x 10 32 (1 - e c ) (l - tcj^A (—) T$e npe cm^s" 1 , (Bl) 

where Tg is the temperature in units of 10 9 K. The step function Q npe is 1 if the direct Urca process is allowed (see 
Ref. [22[), it is otherwise. The factor (A/) is a statistical factor that depends on the phase of nuclear matter under 
consideration. When all particles are in the normal phase, one finds that 

A/ = A/ = — - . (B2) 

60 v ; 

When neutrons or protons are superfluid this quantity is multiplied by a reduction factor, Rx, and one has that 

AI = AI Q R X . (B3) 

Here we consider only nn pairing, and the index X = A, B, C, depending on whether neutrons pair in the 1 Sq channel 
(which corresponds to X = A), in the 3 P2(mj = 0) channel (which corresponds to X = B), or in the 3 P2(mj = 2) 
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channel (which corresponds to X = C). In Ref. |23| it is found that for temperatures smaller than the critical 
temperature of the corresponding supcrfluid phase, the reduction factors are well approximated by 



3..1 



R A = (0.2787 + v/(0. 7213) 2 + (0.1564w A ) 2 ) exp 2.9965 - J(2.9965) 2 + v A , (B4) 



R B = 0.2854 + y/(0. 7146) 2 + (0.1418i; B ) 2 exp 2.0350 - J(2.0350) 2 + u| , (B5) 



0.5+(0.1086w c ) 2 / / ; ttt\ , s 

Rc = l + (0.2347, c )2 + (0 2023, c)4 + °" 5 ^ (' ~ ^^) < ( B6 ) 

where vx = £ T , and Ax (T) is the temperature dependent supcrfluid gap of the phase under consideration. For 
temperatures larger than the critical temperature the reduction factors are equal to one. It is worth noticing that 
while in the phase A and B the reaction rates are suppressed exponentially at low T, the reduction factor for the phase 
C varies as T 2 . The case when both neutrons and protons are in a superfiuid phase is also considered in Ref. |23| . 
where the reduction factors are also numerically computed. 
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